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Time-resolved Faraday rotation measurements have proved transformative in the investigation of 
spin dynamics in semiconductors. In materials with spin lifetimes which are on the order of, or 
greater than, the laser repetition time, the collective effect of spin polarization due to the whole 
pump pulse train becomes important. Here, we discuss a relative phase shift which results from 
these spins. We derive and experimentally validate a closed-form expression which describes this 
phase shift, and characterize it throughout parameter space. A spin lifetime measurement based on 
this phase shift is described, and we discuss situations in which the model used must be augmented 
to be applicable. 
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I. INTRODUCTION 


Time-resolved Faraday rotation (TRFR) measure¬ 
ments have proved to be a powerful tool for the study 
of dynamical spin phenomena in semiconductors^— . 
In its most basic form, TRFR is an optical pump- 
probe technique that can be used to monitor the 
temporal evolution of an optically-injected spin po¬ 
larization and determine the electron g-factor and 
spin lifetime in materials with advantageous selec¬ 
tion rules^. Variations on TRFR measurements, 
including spin-drag and Larmor magnetometry, en¬ 
able the study of spin transport properties^, spin- 
orbit fields^, dynamic nuclear polarization^—, 
and many other phenomena. When the spin polar¬ 
ization lifetime T 2 * is on the order of or larger than 
the laser repetition time, the spins excited by previ¬ 
ous pump pulses make non-negligible contributions 
to the Faraday rotation signal, resulting in resonant 
spin amplification (RSA)^2r—. In samples with long 
lifetimes, RSA can provide a more accurate way to 
measure T 2 * than TRFR—. In this paper, we inves¬ 
tigate another effect of previous pump pulses which 
results in an effective phase shift that depends on 
the applied magnetic field. This phase shift must be 
accounted for in any situation in which the contri¬ 
bution to TRFR from previous pulses is expected to 
be important, and can provide another method for 
determining the spin lifetime. 


II. TRFR MEASUREMENTS 

In a typical TRFR measurement, a circularly po¬ 
larized laser pulse excites a spin polarization parallel 
to the direction of laser propagation, which we take 
to be z. After some time delay, a linearly polarized 
probe pulse interrogates the system by measuring 
Faraday rotation, which is proportional to the spin 
polarization along z. Typically, an ultrafast laser 
is used to generate the pulses and the time delay 
between pulses is controlled with a mechanical de¬ 
lay line. In this way, the shortest accessible time 
scale is set by the pulse time and the precision with 
which the delay line can be positioned. The spin 
polarization lifetime can be found by measuring the 
timescale of the decay of Faraday rotation as At, the 
delay between pump and probe pulses, is swept. For 
a measurement of the g-factor, a magnetic field is 
applied in the Voigt geometry along x. Spins will 
then precess in the y-z plane at their Larmor fre¬ 
quency ULl = gtisB/h, where g is the electron g 
factor, y,B is the Bohr magneton, and h the reduced 
Planck constant. In this geometry, a single pump 
pulse gives rise to a spin polarization of: 

sb = Acos(HiAt) exp[(— At)/T£ ]0(At) (1) 

where A is the initial spin polarization and ©(At) is 
the Heaviside step function which enforces causality 
in the expression. 

When the spin lifetime is such that an apprecia¬ 
ble spin polarization remains after the laser repeti¬ 
tion time, spins due to previous pump pulses must 
be taken into account as well. It is convenient at 


this point to use phasor notation, s —> s where s is a 
complex spin phasor. In this notation, Sz is under¬ 
stood to be the real part of s, Re(s) and Sy becomes 
the imaginary part of s, Im(s). In this case, each 
pump pulse contributes: 

s = Aexp[(ifi L - l/T 2 *)At] 6(At) (2) 

to the total spin. Summing over previous pulses, the 
total s becomes: 

s = A exp [(iO.L — V^ 2 *)(Ai + nT rep )] 

n 

x 0(A t + nT rep ) (3) 

where T rep is the laser repetition time. 


A. Experimental data 

Figure 1 (a) shows experimental data taken on a 
commercially purchased sample consisting of a 2 ym 
layer of n-GaAs doped at 1 x 10 1 ' cm~ 3 atop a semi- 
insulating GaAs substrate. Measurements were per¬ 
formed at 10 K as a function of external magnetic 
field B ext and pump-probe delay time At. A mode- 
locked Ti:Sapphire laser tuned to A = 818.2 nm with 
a pulse time of approximately 3 ps and a repetition 
time T rep of 13.2 ns was used. Pump and probe 
beam powers were measured to be 500 /zW and 75 
/zW, respectively, with beam diameters of approxi¬ 
mately 35 /zm each. Here, the magnitude of the sig¬ 
nal shown is in millivolts as measured on a lock-in 
amplifier, which is proportional to 9f, the Faraday 
rotation angle, but taken to be an arbitrary unit. 
The data set was then fit simultaneously to extract 
the electron g factor and spin lifetime. The data are 
fit to the real part of Eq. (3), with a term added 
which takes into account the contribution from op¬ 
tically excited carriers before recombination, which 
decays rapidly. In total, the fit equation used is: 

9f = Re| y^Aexp[(zfl£ — l/T 2 *)(At + nT rep )\ 

n 

x 0(At + nT rep )} ^ 

+Ao cos(HiAt) exp[At/To] x 0(A t) 

The fit finds g = —0.3788^, T 2 * = 10.2 ns, and op¬ 
tical spin lifetime To = 87 ps. The ratio of Faraday 
rotation from long-lived spins to spins from optical 
carriers, A/Ao, is found to be 2.0. Figure 1 (b) 
shows computed data using the results of the fit. 
Both plates in Fig. 1 have been normalized using 
the same normalization factor. 
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FIG. 1. Normalized lock-in signal (proportional to Fara¬ 
day rotation angle #/) as a function of delay time At 
(horizontal axis) and external magnetic field B ext (ver¬ 
tical axis) for experimental data (a) and data generated 
using Eq. (4) with fit parameters from experimental data 
(b). 


B. Resonant Spin Amplification 


Measurements which rely on resonant spin ampli¬ 
fication take advantage of the collective effect of spin 
polarization resulting from previous pulses. Typi¬ 
cally, the external magnetic field is scanned instead 
of the delay time. When the magnetic field reaches 
a value at which spins precess an integer number 
of times in the laser repetition time, spins result¬ 
ing from all previous pump pulses are aligned along 
z when the next pulse arrives, resulting in a large 
spin polarization at At = 0 + . When a half-integer 
number of precessions occur, spins antialign and a 
smaller polarization is seen. This is shown diagra- 
matically in Fig. 2 plates (a) and (b). These figures 
show the y and z components of spin polarization at 
At = 0+, immediately after a pump pulse arrives, 
normalized to the initial spin polarization of a sin¬ 
gle pulse A. Contributions to the total polarization 
from five additional pump pulses are shown in blue, 
along with the total spin polarization in red. Here, 
data was synthetically generated using g = —0.3788, 
T rep = 13158 ps and T 2 = 2 T rep . 

The degree to which the spin polarization at At = 
0 + exceeds A when a resonance condition is met is 
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FIG. 2. Orientation of total spin polarization (red) and 
contributions due to the five most recent pump pulses 
(blue) for maximum and minimum amplitude [(a) and 
(b), respectively] and minimum and maximum phase [(c) 
and (d), respectively] immediately following the arrival 
of a pump pulse (At = 0 + ). In plates (c) and (d), the 
green arrow indicates the sum of all previous pulses, ex¬ 
cluding the most recent pulse. Data generated using Eq. 
(3), with T£ = 2T rep . All vectors are normalized to 
A = 1, and inset numbers label lines of constant radius. 


determined by the lifetime of the spins; longer life¬ 
times result in higher and sharper resonance peaks 
as measured by scanning the external magnetic field. 
The g-factor and lifetime can be found by fitting this 
data to Eq. 3. This method can give better fits of 
the lifetime than standard TRFR scans when it is 
large compared to the accessible delay times^. 


III. PHASE SHIFT DUE TO PREVIOUS PULSES 


By pulling terms out of the sum in Eq. (3) which 
do not depend on n and enforcing At > 0, s becomes: 

s = Aexp[(iflL — 1/T 2 * )At] 

x ^2 exp[(*OL - l/T 2 *)nT re p] (5) 

n> 0 


where now the sum is in the form of a converg¬ 
ing geometric series with common ratio exp[(if 2 £ — 


FIG. 3. Effective phase <j> (a) and amplitude r (b) as 
defined in Equations (7) and ( 8 ) from model data with 
a range of lifetimes (solid lines) and experimental data 
(black squares). 


l/T 2 )T rep ]. The sum can therefore be evaluated: 


E G 0 -)" = <6) 

n> 0 


with 0 = f2iT rep and x = T rep /T £. This can be 
rewritten as a complex number re^ whose ampli¬ 
tude r and phase are functions of the external 
magnetic field and spin lifetime: 


r = (l — 2e~ x cos( 0 ) + e~ 2x ) * 


4> = tan 1 


/ e~ x sin( 0 ) \ 
\ 1 — e~ x cos( 0 ) / 


( 7 ) 

( 8 ) 


By taking the real part of Eq. (5) using these ex¬ 
pressions, Sz becomes: 


Sz = rAcos^iAt + (j))e At / T 2 ( 9 ) 


These expressions are functionally equivalent to the 
closed-form expression found in Ref. 17, but are 
valid for At € [0,T rep ) instead of At € [— T rep , 0), 
and are written in a form that makes explicit the fact 
that previous pulses contribute an effective phase 
shift <fr to TRFR measurements in addition to mod¬ 
ifying the amplitude of the spin polarization. This 
phase shift is equivalent to tan^ 1 (sp/s£)| z ^ t=0 +. 

In an effect analogous to peaks in RSA amplitude, 
there exist values of the external field where the sum 
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FIG. 4. Path traced out by total spin polarization 
s(At = 0 + ) as applied magnetic field is swept. In the 
limit T£ —y oo, the radius of the circle generated goes to 
infinity, which results in maximum and minimum phases 
of 7 t/ 2 and — 7r/2, respectively. 


of the spin polarization due to previous pulses will 
result in a maximum effective phase. At these field 
values, the spin polarization immediately preceding 
the arrival of a pump pulse s( 0 _ ), will lay in the 
sample plane along This condition is shown in 
Fig. 2 plates (c) and (d). In these figures, the green 
vector indicates the sum of previous pulses exclud¬ 
ing the most recent pulse. Again, the red and blue 
arrows show the total spin polarization and the first 
five terms in the sum, respectively. At field values 
which satisfy this condition, the effective phase </> in 
Eq. (9) is maximized. 

Figure 3 shows plots of r and (j> versus external 
magnetic field calculated using Equations (7) and 
( 8 ) with a range of lifetimes. Individual delay scans 
taken from the experimental data in Fig. 1 (a) were 
fit to Eq. (9) to extract r and 0, and the results are 
plotted as black squares. Calculated values of r and 
cf> using the experimentally determined lifetime and 
g factor from the simultaneous fit above are shown 
as a black line to highlight the close agreement with 
the model. As the spin lifetime increases, the max¬ 
imum and minimum observed phase angles increase 
in magnitude and are found closer to the external 
magnetic field values corresponding to an RSA peak. 


As the magnetic field is swept, the total spin po¬ 
larization at At = 0+, depicted by the red arrows in 
Fig. 2, traces out a circular path in the y-z plane. 
This path is shown in Fig. 4 for the lifetimes used 
previously in Fig. 3. As the lifetime grows, the 
radius of the circle increases, and in the limit that 
T 2 * goes to infinity the radius of this circle diverges. 
Accordingly, the maximum and minimum phase an¬ 
gles approach 7 t/2 and — 7 r/ 2 , respectively. The rate 
at which s(Af = 0 + ) moves around the circle as the 
magnetic field is swept is not constant, but rather ex¬ 
hibits a rapid shift in phase as the field value moves 
past an RSA peak and then a slow shift in phase 
through the RSA minimum. 

A. Measurements using <f> 

Since arctangent grows monotonically with its ar¬ 
gument, the extrema of (j) can be found by extrem- 
izing the argument of the arctangent in Eq. ( 8 ). 
Setting its derivative equal to zero and solving gives 
the values of 0 which extremize <fi: 

®extr = 2 wr±cos 1 (e x ) ( 10 ) 

As T 2 grows (causing x to shrink with constant 
T rep ), the values of @ extr grow closer to 2 n 7 r, which 
correspond to peaks in r. Substituting Eq. (10) with 
n = 0 into Eq. ( 8 ) gives the maximum phase angle 
as a function of x, given by: 



which is plotted in Fig. 5 (a). 4> m ax has a monotonic 
dependence on the spin lifetime, and can therefore 
be used to measure T 2 when it is sufficiently large 
compared to the laser repetition time. Such a mea¬ 
surement could be performed by first setting At to a 
small negative value and tuning the magnetic field so 
that the Faraday rotation disappears. At this mag¬ 
netic field, the polarization due to previous pulses 
will lay in the sample plane, and correspond to a 
condition of maximum or minimum phase. Then, a 
delay scan would be performed and fit to extract (j), 
which, by inverting Eq. (11), can be used to find T|. 
In the experimental data shown in Fig. 3, the maxi¬ 
mum observed phase angle was 16.4°, which gives a 
lifetime of 10.4 ns, close to the value obtained above. 

TRFR can be used to measure the local magnetic 
field when the g factor is known by a process dubbed 
Larmor magnetometryiS. Traditionally, a delay scan 
is performed over at least a few periods of preces¬ 
sion and a fit to the precession frequency extracts 
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FIG. 5. (a) Maximum phase angle (j>max versus T£/T rep , 
calculated using Eq. (10). (b) Delay time corresponding 
to the 15 th zero crossing of Ss. (vertical axis) for a subset 
of the lifetimes defined in the Fig. 3 legend versus exter¬ 
nal magnetic field. (Inset) Zero crossing delay time for 
the 6 th crossing. In this range of external magnetic field 
strength, a breakdown of the one-to-one correspondence 
between Ato and B occurs for long spin lifetimes. 

the magnetic field present. The speed with which 
such a measurement can be performed is set by the 
time it takes to complete the delay scan, typically on 
the order of 10 to 100 seconds. Measurements have 
also been performed by maintaining a fixed delay 
and continuously monitoring Op — . This technique 
does not generally allow for a good determination of 
the magnetic field because Op is not single-valued 
when the magnetic field changes, and when near 
a local extrema, Op changes slowly with magnetic 
field, reducing or eliminating the ability to measure 
changes in the magnetic field. Additionally, shifts 
in laser power or slow drifts in pump-probe overlap 
are indistinguishable from changes in magnetic field. 
These limitations can be avoided if instead the delay 
time is continuously adjusted to follow a zero cross¬ 
ing of Op of known order. The position at which 
a zero crossing occurs is unaffected by changes in 
laser power or pump-probe overlap, and, with care, 
a zero crossing can be picked whose delay time has 
a one-to-one correspondence with the magnetic field 
strength over a range of values. A measurement can 
then be completed in the time it takes to adjust the 
position of the mechanical delay line once and take 
a measurement of Op, which can take place in ap¬ 


proximately one second. 

A zero crossing will occur when s* = 0, or from 
Eq. (9), the n th crossing occurs when fl^At + <f> = 
(2 n — l)ir/2. The presence of 4> in this expression 
complicates this measurement. Figure 5 (b) shows 
the delay time at which the n = 15 crossing occurs 
versus magnetic field for a range of lifetimes. Again, 
colors refer to the legend in Fig. 3. When the spin 
lifetime is long, the change in phase due to previous 
pulses is large, and results in a series of sharp steps 
in the delay time at which a zero crossing occurs 
that correspond with RSA peaks as the magnetic 
field strength is scanned. It can be shown that if 
a zero crossing occurs at a delay time that satisfies 
At > T rep / 2, there will always be a one-to-one re¬ 
lationship between the local magnetic field strength 
and the delay time at which the crossing occurs. The 
inset in Fig. 5 (b) shows what happens when this 
relation is violated for a subset of the spin lifetimes 
defined in Fig. 3. Here, the 6 th crossing has been 
chosen. In this situation, for long spin lifetimes, the 
local magnetic field can no longer be uniquely de¬ 
termined from the delay time of the zero crossing. 
In order to measure the local magnetic field using 
this technique, the spin lifetime must now be known, 
which adds uncertainty to the measurement of the 
magnetic field. Additionally, it must be known ex¬ 
actly which crossing is being followed. This could 
be ascertained by performing a full delay scan at a 
fixed, known magnetic field. 

B. Breakdown of the model 

In the preceding analysis, it is assumed that the 
contributions from previous pulses decay exponen¬ 
tially with the time elapsed since the pump pulse 
arrived at the sample. In situations where there is 
electron drift or diffusion, each term in the sum in 
Eq. (3) must now include a factor which takes this 
into account. In general, this will depend on many 
parameters, including but not limited to beam size 
and shape, mobility, in-plane electric fields, and dif¬ 
fusion constant. Equation (5) becomes: 

s = £ {Aexp[(iflz, - l/T 2 *)(Af + nT rep )] 

n 

x F n (At)Q(At + nT rep )} (12) 

where F n (At) takes into account changes in overlap. 
Since F n (At) depends on n and At, it cannot be 
factored out of the sum. As a result, this expres¬ 
sion cannot be simplified to the form of Eq. (9), 
and the results of a fit to that functional form will 
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no longer directly give the precession frequency or 
lifetime. Care must be taken to characterize and 
correct for this effect whenever the in-plane motion 
of electrons is expected to be important. Also, we 
assume that spins resulting from each pump pulse 
are not affected by subsequent pulses. If the optical 
power is sufficient that the number of optical car¬ 
riers is on the order of or greater than the number 
of dopant carriers within the beam area, this is no 
longer a good assumption^, and a factor which ac¬ 
counts for this must be included in the sum. This 
factor will not depend on At, however, so that Eq. 
(9) may still be valid, though the expressions for r 
and 4> must be revised. This effect can be minimized 
by decreasing the pump and probe powers. 


IV. CONCLUSION 

In conclusion, we have recast expressions used in 
the context of RSA measurements in a form that 
highlights an effective phase shift due to the collec¬ 
tive action of spin polarization from previous pulses 
in TRFR measurements. We experimentally mea¬ 
sure this phase shift in n-GaAs and find that it is in 
good agreement with the model. In situations with 
long spin lifetimes, the phase shift can become a sig¬ 
nificant factor in TRFR measurements. We suggest 
a measurement which utilizes the phase to provide 
another means of accessing the spin lifetime when it 
is long compared to the laser repetition time in which 
At is scanned instead of the external magnetic field. 
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